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GYRO-SYNCHROTRON EMISSION AND 
ABSORPTION IN A MAGNETOACTIVE PLASMA 

Reuven Ramaty* 

ABSTRACT 

The intensity, spectrum and polarization of gyro- synchrotron radiation from 
an arbitrary distribution of electrons in a magnetoactive plasma are calculated. 
These calculations are based on detailed theories of gyro- synchrotron emission 
and absorption as well as on the equations of radiation transfer in the limit of 
large Faraday rotation. By performing extensive numerical calculations, for 
both isotropic and anisotropic electron distributions, we investigate the influence 
of reabsorption and the Razin effect on the spectrum and polarization of the ra- 
diation. Using the same formalism and calculations, we also investigate the 
possibility of negative absorption in an anisotropic gyro- synchrotron source. 
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GYRO-SYNCHROTRON EMISSION AND 
ABSORPTION IN A MAGNETOACTIVE PLASMA 

I. INTRODUCTION 

Gyro- synchrotron radiation is the electromagnetic emission generated by 
mildly relativistic electrons moving in a magnetic field, while the term synchro- 
tron radiation is used to describe the emission from ultrarelativistic electrons. 
The theory of synchrotron emission and absorption, in vacuum as well as in a 
plasma, has been treated by a number of authors (e.g. Ginzburg and Syrovatskii, 
1969), and it is widely believed that it can account for the observed radio emission 
from a variety of cosmic sources. Synchrotron theory, however, is not directly 
applicable to radio emission produced by mildly relativistic electrons in regions 
of strong magnetic fields and dense plasmas. In such regions, the ambient 
plasma and gyro frequencies are large and comparable, within an order of mag- 
nitude, to the observed radio frequencies; the propagation of electromagnetic 
waves is anisotropic; there is more than one distinct mode of wave propaga- 
tion; and there may be a large Faraday rotation. This type of radiation is prob- 
ably responsible for radio emission from large solar flares as well as for emis- 
sion generated in the atmospheres of certain stars and planets. For the detailed 
understanding of its spectrum and polarization, it is necessary to consider in 
detail the theory of gyro- synchrotron emission and absorption in a magneto- 


active plasma. 


The theory of gyro-synchrotron emission, for electrons in circular orbits 
in vacuum, was first treated by Schott (1912). Subsequently, Schwinger (1949), 
Takakura (1960a) and Landau and Lifshitz (1962) have further developed this 
theory, and Takakura (1960b) has extended it to electrons in helical orbits. The 
general problem of electromagnetic radiation by charged particles in a mag- 
netoplasma has been treated by a number of authors (Twiss, 1958; Eidman, 1958; 
Liemohn, 1965 and Mansfield, 1967). These treatments provide expressions for 
the spectral and angular distribution of the radiation, in a given mode of wave 
propagation, for single electrons of given energies and pitch-angles. By con- 
sidering the special case oi in isotropic ambient medium, Ramaty (1968) ob- 
tained an expression for the total radiation spectrum, integrated over all emis- 
sion angles, for electrons in circular orbits. A numerical analysis of this 
spectrum revealed that if the Lorentz factor of the radiating electron is less 
than the plasma -to-gyro-frequency ratio of the ambient medium, the radiation at 
low frequencies is strongly suppressed. This effect is well known for synchrotron 
radiation, and is generally referred to as the Razin effect. 

The radiation at low frequencies, however, is influenced not only by the 
Razin effect, but also by a variety of propagation phenomena, such as absorption 
below the plasma frequency, gyro- resonance and free- free absorption by the 
ambient electrons, and gyro- synchrotron reabsorption by the radiating electrons 
themselves. 
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The theories of gyro- resonance and free-free absorption were discussed 
by several authors (e.g. Kundu, 1965), and we shall not treat them, therefore, 
in the present paper. The theory of gyro- synchrotron reabsorption was treated by 
Twiss (1958), for radiation in vacuum, and was extended by Kawabata (1964) to 
radiation in a dilute plasma with Faraday rotation. An attempt has been made 
by Kai (1965) to apply Kawabata's (1964) results to solar radio emission, but 
because of the complexity of the equations, no detailed numerical results were 
obtained. In a complete treatment of gyro-synchrotron emission and absorption, 
the Razin effect and reabsorption must be treated in a consistent fashion since 
these two processes are highly interdependent. As the influence of the ionized 
medium becomes more pronounced, the Razin effect not only suppresses the 
emissivity of individual electrons but also reduces the corresponding absorption 
coefficient. Such a simultaneous treatment of these two effects is not available 
in the published literature. 

Closely connected with the theory of gyro- synchrotron reabsorption is the 
problem of negative absorption. The possibility of the amplification of synchrotron 
radiation was first considered by Twiss (1958), but Bekefi, Hirshfield and Brown 
(1961) showed that under the circumstances described by Twiss (1958) (isotropic, 
ultrarelativistic electrons in a vacuum) the absorption coefficient will always be 
positive. Negative absorption, however, may occur if the radio source consists 
of mildly relativistic electrons. Bekefi, Hirshfield and Brown (1961) have in- 
vestigated this possibility and showed that if the electron spectrum increases with 
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increasing energy but cuts off before the electrons become ultrarelati viatic, 
maser action will occur just below the gyro-frequency and its second harmonic. 
Even though Twiss' original suggestion for maser action was proven to be unten- 
able, the possibility of amplification of synchrotron radiation from ultrarelativistic 
electrons was reconsidered recently by a number of authors. McCray (1966) and 
Zheleznyakov (1967) have shown that maser action may occur if the radiation is 
produced in a plasma and the electrons have an inverted energy spectrum. 
Heyvaerts (1968) has considered the possibility of maser action due to pitch- 
angle anisotropies and showed that amplification of synchrotron radiation is 
possible if the angular spread of the particle distribution is of the order of the 
aperture angle of the radiation pattern. Obviously, for ultrarelativistic electrons, 
extremely well collimated particle beams are necessary. For mildly relativistic 
electrons, however, maser action would require much smaller anisotropies. The 
possibility of the amplification of gyro-synchrotron radiation owing to an anistropic 
pitch-angle distribution, has not been yet considered in the published literature. 

In the present paper we shall provide a detailed treatment of the theory of 
gyro-synchrotron emission and absorption by an arbitrary distribution of radiat- 
ing electrons in a cold, collisionless magnetoactive plasma. We shall evaluate 
Stokes' parameters by using a general expression for gyro-synchrotron emis- 
sion (Liemohn, 1965), a detailed formula for the absorption coefficient (Bekefi, 
1966), and the theory of gyro- synchrotron radiation transfer in the limit of large 
Faraday rotation. By performing detailed numerical calculations we investigate 
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the simultaneous effects of the ambient medium and reabsorption. We find that 
both these effects tend to suppress the radio flux at low frequencies but lead to 
different polarization of the resultant radiation. We also consider the effects of 
an anisotropic electron distribution. Since the radiation pattern of ultrarelativistic 
electrons is much narrower than that of mildly relativistic particles, the emis- 
sion spectrum and polarization from an anistropic distribution will strongly de- 
pend on the direction of observation. We find that if this direction differs from 
that of maximum anistropy, the radiation at high frequencies is strongly sup- 
pressed while the emission at the longer wavelengths remains essentially un- 
changed. This suppression effect coupled with the Razin effect at low frequencies 
may lead to the total suppression of the entire radiation spectrum. 

Finally, we consider the possibility of maser action in an anistropic gyro- 
synchrotron source. We find that maser action is indeed possible, even if the 
electron spectrum decreases with increasing energy, and, for moderate anistropies 
it occurs at essentially the same frequencies as found by Bekefi, Hirshfield and 
Brown (1961) for the isotropic distribution with an inverted energy spectrum. 

E. GYRO-SYNCHROTRON EMISSION AND ABSORPTION 

We consider a homogeneous system of energetic electrons moving in a 
homogeneous, cold and collisionless electron plasma, permeated by a static 
and uniform magnetic field. If the density of the plasma electrons is much higher 
than that of the energetic particles, wave propagation in this system is governed 
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solely by the ambient medium. Waves of two modes can propagate In such a 
medium, the ordinary and extraordinary modes of magnetoionic theory (Ratcliffe, 
1959; Ginzburg, 1964). The indices of refraction, n ± , and polarization coefficients, 
a^ ± and a^ ± , as functions of the frequency v, and of the angle 8 between the wave 
normal k and the static field B, are given by (Ginzburg, 1964): 


n t 0) ~ 1* 






± fv 4 v* sin* t? ( l ' p 2-I ' i ) 2co8i - v 3 ) -v a sin a 6 

2v(v p 2 - v 2 ) cos 8 


( 1 ) 




v B sin 2 8 ± |i/ 4 v B 2 sin 4 8 + 4v 2 (v p 2 - v 2 ) 2 cos 2 0^* 


/ 2 


( 2 ) 


a^(v,8) = 


J'p 2 v B v sin 8 - a f)t v o 2 i^ B 2 cos 8 sin 8 


e t p b 


V p 2 ^V B 2 COS 2 0 * v 2 ) “ l' 2 [l/g 2 - V 2 J 


(3) 


The subscripts (+) or (-) correspond to the ordinary or extraordinary modes, 
respectively, and the plasma and gyro frequencies, v p and v B , are given by 

_ / / « i/2 _ 1 bB 

% - e(n e /7rm) 1/2 ; " 17 me 


where n e is the electron density and Bis the magnitude of the static magnetic 
field. The polarization coefficients are defined in terms of the components of 
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the electric vector of the radiation field 


h . h. 

ia 0 - r ; iQ k - 

— ■+ 

where E x and E e are the transverse components of E, with E d lying in the plane 
defined by Band k, and E k is the component of E along the wave normal direc- 
tion fc. 

The indices of refraction, n ± , have several cutoffs (n = 0) and resonances 
(n -<»), and for a non-vanishing B there are frequency bands below the plasma 
frequency in which either or both the ordinary or extraordinary mode indices 
are real and greater than unity. It can be shown, however, (Pawsey and Brace- 
well, 1955), that for 9^0, radiation at a frequency v, in the ordinary or extra- 
ordinary modes, will not escape from a radio source for which both n e and B 
fall off toward the observer, unless v is greater than v or v xt respectively, 
where v x ~ (> p 2 + v z/^) 1/ * + V B / 2 - Furthermore, for v > v pt n + is real and 
less than unity, and for v > v %t n_ is also real and less than unity. We shall limit 
our discussion, therefore, to these frequency ranges. 

As can be seen from Equations (3) and (4), the electric vector of the radia- 
tion field in a plasma has in general a non-vanishing longitudinal component, E k . 
Unless we invoke coupling between longitudinal oscillations and transverse 
electromagnetic waves, the energy associated with this component will not escape 
from the radio source, and therefore, in the subsequent calculations we have set 
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a k = 0. We wish to point out, however, that this cannot introduce any serious 
error in our calculations, since from the numerical evaluations which will be 
described below, we find that for the frequency range of interest (v > v p and 
v > v x ), the emissivity associated withE k is negligible in comparison with the 
total radiated energy. 

The polarization coefficient, a ±t for an arbitrary direction of propagation, 

d 

is finite and greater than zero. The radiation is therefore elliptically polarized. 
For the special cases of longitudinal (9 = 0) or transverse (9 -v/ 2) propagation, 
the radiation becomes circularly or linearly polarized, respectively. In partic- 
ular, for 9 - 77/2 

n + 2 = 1 - 

n_ = 1 + 

and 


We now consider the radiation produced by the energetic electrons. If we 
set a k = 0 , the frequency and angular distribution of the emission from a single 
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electron, into the ordinary or extraordinary modes, is given by (Liemohn, 1965): 



where 

sn ± /Ssin </>sin 8 

x s - 1 - n t /Scos <£cos 8 ' (8) 

In these expressions, p is the electron velocity in units of c; y - (l - / 3 2 ) - 1/2 ; 

4> is the electron pitch- angle with respect to B and J s is a Bessel function or 
order s. Equation (7) can be used for all frequencies for which n 2 > 0. Positive 
values of s correspond to the normal Doppler effect, negative values of s, to the 
anomalous Doppler effect and s = 0, to Cerenkov radiation. For values of v 
greater than v p and v x , however, 0 < n 2 <1, and, therefore, the index s ranges 
only over values greater than zero. 

The emissivity v, given in Equation (7) is equal to the energy lost by the 
particle due to radiation into unit frequency interval and unit solid angle and not 
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to the observed radiation by a distant observer. It was shown by a number of 
authors (Epstein and Feldman, 1967; Scheuer, 1968; Ginzburg and Syrovatskii, 
1969) that if the radiating electron moves in a helical orbit, these two expres- 
sions are not identical, the total energy loss being smaller by approximately a 
factor of sin 2 4> than the integrated radiation crossing a large but fixed surface 
surrounding the moving charge. This difference is the result of the nonstationary 
nature of the radiation field which leads to a change in field energy within a fixed 
volume of space. It was shown by Scheuer (1968) and by Ginzburg and Syrovatskii 
(1968), however, that in order to compute the radiation from a given number of 
electrons in a fixed volume in space, one still has to use Equation (7) for the 
power emitted by an electron and not the modified expressions that would cor- 
respond to the observed radiation by distant observers. Equation (7) should be 
treated, therefore, as the volume emissivity from a monoenergetic distribution 
of electrons with unique pitch-angles, normalized to 1 electron of Lorentz factor 
7 and pitch- angle <f>. 

The emissivities, j ± , and absorption coefficients, K ± , for an arbitrary 
energy and pitch-angle distribution can be obtained from the following general 
expressions (Bekefi, 1966; Melrose, 1968): 

j t (v, 6) ~ f V± (v, 0, 7 , <p) f(p')d 3 p' ergs/(sec sr Hz cm 3 ) (9) 
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K ± (v, 6 ) 


( 10 ) 


~ ^7 | V ± ( V< 6 ' y ' TuT [ f (P) “ f (p')] d3 p' C" 1 ' 1 


where f(p)d 3 pis the number of electrons per unit volume, with momenta in 
d 3 p around p, and the sign of K is defined such that by emitting a photon of en- 
ergy an electron changes its momentum from p' to p. Since the photon energy 
hv is much smaller than the electron energy 7 me 2 , f(p) - f(p') can be reduced 
to 


_ _ hv df <9 f 

f(p) - f(p ) - - Pyjj, + tg (ii) 

where A<£ is the change in the pitch-angle of the radiating electron resulting from 
the emission of a photon of energy hv into the direction 6, and is given by 

. hv nfi cos 6 - cos 4> 

tgA<£ - p C /3sin7 (^) 

As can be seen from Equations (10) and (11), for an isotropic distribution (<9f/ 
d<t> = 0), a necessary condition for negative absorption isdf/dy > 0. For an 
anisotropic distribution, however, negative absorption can occur even if <?f /dy < 0. 
A necessary condition for this is 




0 


11 


which simply means that owing to the emission of a photon, the electron pitch- 
angle is changed from a more densely populated to a less densely populated 
pitch-angle state. We shall further investigate this type of instability in the 
numerical calculations to be described below. We wish to point out, however, 
that maser action from a moderately anisotropic distribution will be confined 
mainly to the low harmonics of gyro-synchrotron emission, since for ultra- 
relativistic electrons the emission is strongly peaked in the direction of the 
instantaneous velocity, and amplification will occur only if the angular spread 
of f is of the order of the aperture angle of the radiation pattern. This point 
was made by Heyvaerts (1968). 

We now consider a radio source which has a total volume V and contains N 
electrons with energies greater than some given minimal value. We assume 
that the distribution function f(p) can be separated into an energy dependent 
part, u(y) , and a pitch-angle dependent part, g(4>) such that 

N 1 dy 

f(P) = \TI7 d^ u OO g ( 0 > • ( 13 > 

In terms of these functions, j and K can be written as: 


j ± (*, e ) 


2tt 


r r 1 

J u(y) dy J 


d cos <p g{4>) 7} (v, 6, y, 4> ) 


ergs 

cm 3 sec sr Hz 


(14) 
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K± (v, 0) 


2tt 


mv 2 n 


r r 1 

I u(y)dy I dcos <£g(<£) rj ± (v , 6, y, <p) 

J x J- 1 


[iy 2 d /u (y)\ n/3 cos 8 - cos 4> 1 dg(0) 

u ( y) \j3y 2 ) + y{3 2 sin^ 8W 


cm' 1 . (15) 


By substituting Equation (7) for v± we can directly evaluate j and K for both the 
ordinary and extraordinary modes. The double integrals, dyd cos 4>, in Equations 
(14) and (15) can be reduced to single integrals by virtue of the delta function in 
Equation (7). If 6 / tt/ 2 , this eliminates the integral over d cos <fi, while if 
6 = tt/2 the delta function is used to eliminate the integral over dy. Carrying 
out these integrations, we obtain: 



BN 

e 3 / v \ 

(16) 

j± (*. 0) 

" V 


K ± (v. 6) 

N 

0 ( V \ 


‘ BV 

(2T7) 2 eHj— , 9\ . 

(17) 


For 8 / 77 / 2 , the functions G and H are given by 



277 

cos 6 



s gK) / 

1 ~n ± /3cos 6 cos <p s i+ a 2 \ 

e ± 


d - in n ± 
+ 6 v 



sin Ms' ( X s) 


+ 


(cot 0 


-fi 


COS <p\ 

sin 6 ) J,( x s) 



1 

/3y 2 d / u(y)\ n/3 cos 8 - cos 

liCy) dy ^ J y/3 2 sin</> 


1 dg\ 

g (4>) d<fc\ 


( 18 ) 


13 


where 


sn t /3 sin 9 sin 4> B 

1 - n ± yScos 9 cos 4> % (^®) 


Sj 2 - 7(l ± n ± /3 cos #) (20) 

B 

1 - sVg/Tv 

cos4> a - n ± /3cos9 ’ (^ 1 ) 


For 9 = 77/2 we obtain: 


G ± 

« ± 


2tt n ± 


L s ^ 2 

s= 1 



I " 1 

+ T&W J t d cos 7 Y ± g(<£ s ) 


1 

1 ( /v b\ 2 /_ ^» 2 _d_ / u \ 

n ± 2 \ v ) WL 


cot c j> 1 dg 

y fi * g(<£) d0 


) 


( 22 ) 


where 



x § - sfi t m ± sin <t> . 


( 23 ) 


( 24 ) 

( 25 ) 
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In previous studies of the effect of the ionized medium on solar radio emis- 
sion (Ramaty and Lingenfelter 1967, 1968), it was convenient to introduce a param- 
eter a, defined as 

3 V B 

a = JIT ’ (26) 

p 

Then, in terms of this parameter, gyro- synchrotron emission from an electron 
of Lorentz factor y is strongly suppressed at the low frequencies if ay < 1, but 
remains unaffected by the ionized medium if ay £ 1 (Ramaty, 1968). Moreover, 

I if ay << 1, the total emission will also be strongly suppressed. Since the indices 

| of refraction and polarization coefficients, as functions of v / v B , depend only on 

i the "intrinsic" parameters a and 9, the functions G and H also depend on these 

I parameters only and on the normalized distribution functions u(y) and g (<£) , 

but are independent of such "extrinsic" parameters as the intensity of the mag- 
I netic field, the size and shape of the radio source and the densities of the en- 

| ergetic and ambient electrons. 

m. POLARIZATION AND RADIATION TRANSFER 

Partially polarized radiation at a given point in space is completely deter- 
mined by 4 independent parameters such as the Stokes parameters I , Q, U, and 
V. These functions in turn, can be obtained from any other set of Independent 
parameters that specify the radiation. In a magnetoactive plasroa, the natural 
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set of parameters consists of the intensities of the ordinary and extraordinary 


modes, I + and I_, and the phase relations, I t and I c , which are defined in terms 
of I* and the phase difference, 8, between the modes: 

I f = (l + I_) 1/2 sin 8 

I c = (l + l.) V2 cos 8 . 


Since, as mentioned above, the longitudinal component of E is set to zero (a k = 0), 
the radiation field is completely determined by its transverse components. Using 
Equation (4), E x and E g , as functions of I ± and 8, can be written as: 



(27) 


where the unit vectors e ± are defined by 



and where the x and 8 axes are chosen such that the 9 direction is along the 
perpendicular component, E^, of the static field B. In terms of the transverse 
components of the electric field, Stokes' parameters are given by (Born and 
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Wolf, 1964) 


q = £(e„e;-e 9 e;) 

0 = (*.*,’ **>*.") 

v = et* (e,e;-e 9 e;) 


•\ 


> 



( 28 ) 


The degree of polarization, v, the ratio of the axes of the polarization ellipse, p, 
and the angle x between the major axis and the x direction can then be written as 


77 = (q 2 +u 2 +V 2 ) 1/2 /l 

p = tg" 1 /?; sin 2/3 = v/(0 2 +U 2 + V 2 ) 1/2 L (29) 


t g 2 x = U/Q 


These quantities are defined such that, if seen by an observer, the electric vector 
rotates clockwise or counterclockwise if V ip negative or positive, respectively. 

V = 0 corresponds to linear polarization and Q = U ■ 0 corresponds to circular 
polarization. The angle x(0 < x <7T )> as seen by an observer, is measured counter- 
clockwise from <;he x direction; ifU>0, 0<x < w /2 and if U < 0, -n/2 <\<-rri it 
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U = 0, * = 0 or * = 7 t/ 2 for Q > 0 or Q < 0, respectively. These definitions are 
consistent with Figure (1.8) in Born and Wolf (1964). 

Substituting Equation (27) into Equations (28), we get 


I = I + + I- ] 



For radiation in vacuum, a ± = ±1, and Equations (30) reduce to 

9 


I = I* + I. 
Q = 2I C 
U = 21, 

V = I + - I. 




> 


J 


(30) 


(31) 


Equations (30) or (31) provide the transformation between the two independent 
set of parameters, I , Q, U, V and I t , I,, I c . Similar equations were obtained by 
Kai (1965). They contained, however, some numerical errors, and therefore 
could not be used in the present study. 
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If the radiation is completely polarized, the phase difference 8 has a well 


defined value, and 


I s 2 + I c 2 = I + I. • 

Using Equations (29) and (30) or (31), we verify that in this case the degree of 
polarization, -n, indeed equals unity. 

For a single particle in vacuum, S = 0, and, therefore, the radiation speci- 
fied by Equations (31) is elliptically polarized with the major axis perpendicular 
to B x . For ultrarelativistic particles at high frequencies I + = I_, and the ra- 
diation becomes linearly polarized with the plane of polarization perpendicular 
to B r For a single particle in a plasma, 8 is not equal to zero and its value 
depends in a complicated way on the energy and pitch-angle of the particle, as 
well as on the direction of propagation and the properties of the plasma. We 
shall circumvent these complications, however by limiting our treatment to the 
case of large Faraday rotation in the source region, for which the phase relations 
are randomized and I s = I c = 0. The assumption of large Faraday rotation in 
the source is almost always well satisfied for gyro-synchrotron radiation. The 
rotation angle, A </>, is approximately given by 


where L is of the order of the linear dimensions of the system and K is the wave- 
length of the radiation. Since gyro-synchrotron radiation generally occurs at 
frequencies v, comparable, within an order of magnitude, to both v p and , 
whereas L » k. Ay; will in general be much larger than 2v and the amount of 
Faraday rotation will be large. 

In order to obtain the intensity and polarization of the radiation from a 
system of particles, it is necessary to solve the transfer equations for either 
I,Q, U and V or I ± , I s and I c . General treatments of such equations, without 
specific solutions, were given by Kawabata (1964), Kai (1965), and Ginzburg and 
Syrovatskii (1969). The phase relations, I s and I c , provide coupling between 
the two modes, and this considerably complicates the solution of these equations. 
If the Faraday rotation is large, however, I g = I c = 0, the transfer equa- 
tions for I + and I_ decouple and take on the simple form 

dl + (v, 6 ) 

~ "T z + K ± (v, 6) I ± (v, 6) j ± 0, 0) (32) 

where j ± and K ± are given by Equations (9) and (10). 

For a homogeneous source region, the solution of Equation (32) is given by 

I ± (v,6) = jr^ [l -exp (-K ± L)] (33) 
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where L is the depth of the source. If the radio source has an effective area A 

such that its volume is V * AL, and if in addition to the suppression of the 

longitudinal component of the radiation field, discussed above, no other bound- 
ary conditions are required, then by using Equations (16), (17), and (33), the 
intensity of the escaping radiation can be written as 

e 2 (BN/A) G ± (vA b ) f / n \1 

I±(V,5) = ('^'mc 2 (N/BA) H, (vA b ) [l ~ exp ^ (2 tt) 2 e BA H ±)J * < 34 > 


For I = I = 0 Equations (30) also reduce to 


I + + I_ 


2 


1 - a 


I + 


1 + a 2 . 


+ I. 


1 - a 2 _ 

e_ 

1 +a 2 

e~ 


U = 0 


I + a . I_ a _\ 

v = 2 '7T^T + r^ 

e * e 


(35) 


Equations (34) and (35), with G ± and H ± given by Equations (16) through (25), and 
a ± given by Equation (2) completely determine the intensity and polarization of 

d 

the radiation. Since G(v/v B ) and H(v/V B ) depend only on the "intrinsic" param- 
eters a and 6 and the normalized distributions u(y) and g(<£) , the "extrinsic" 
properties of the source depend only on N/BA and BN. In particular, the 


21 


reabsorption depends only on the parameter N/BA (l.e. on the total number of 
energetic electrons along the line of sight per unit area and per unit B) . 

We consider now the Stokes parameters given by Equation (35) . Using 
Equations (29), we find that the degree of polarization is given by 

Thus, when both modes escape from the source region, because of Faraday 

rotation, n is always less than unity. Since U = 0 (which is also a result of 

Faraday rotation) one of the axes of the polarization ellipse will be parallel to 

Bj. From the numerical calculations which are described below, we find that the 

values of Q are such that the major axis of the ellipse is parallel or perpendicular 

toB , depending on whether I + is greater or smaller than I_ , respectively. 

Linear polarization will occur if Q - n/2, since then, by Equations (6) and (35) 

V = 0 and Q = I_ - I + . Thus, depending on whether I_ is greater or smaller 

— » 

than I + , the plane of polarization is perpendicular or parallel to B I , respectively. 

In the next section we shall use the formalism developed above to study the 
spectral and polarization properties of typical gyro- synchrotron sources. 

IV. NUMERICAL RESULTS 

We first consider monoenergetic electrons with an isotropic pitch-angle 
distribution, g(<£) = l/4n . Using Equation (16), we have e/aluated the emis- 
sivity j , equal to j + + j_ , for a variety of values of a ,6 and y. The resultant 
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spectra, as functions of v/v e , where 


v c - v B y 2 (36) 

are shown as solid lines in Figure 1 for: a = 0.25; 6 - 70° and 45°; and y = 2, 
7 = 5 and y - 00 . The dashed lines in Figure 1 represent the emissivity j UR , 
evaluated from the formula 


^ UR ( 


, ' l9) = ^^ sin0 (7)(^)( /v 


(37) 


where 



(38) 


and 


7 

9 y 2 - 1 7 2 (39) 

4 - ) 2 J 

Equation (37) appears to be the generalization of the well known emissivity 
formulas foi* ultrarelativistic electrons in helical orbits in vacuum (Ginzburg 
and Syrovatskii, 1964, 1965) and for ultrarelativistic electrons in circular orbits 
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1 


i 


in a medium (Ginzburg and Syrovatskii, 1964, 1965; Ramaty and Lingenfelter, 
1967) As can be seen from Figure 1, for both 1 9 = 70° and 9 = 45°, there is a 
significant discrepancy between j UR and the exact spectra, j , obtained from 
Equation (16). The ratios, j UR y/j , of the spectra shown in Figure 1, can be 
compared with the ratios F ur /f that were evaluated by Ramaty (1968) for the 
integrated emissivity (over all angles 0) from mildly relativistic electrons in 
circular orbits in an ionized medium. From this comparison (see Figure 3, 
Ramaty (1968)) we find that the ratios j UR /j , for 9 = 70®, are almost the same 
as the corresponding ratios F UR ^F . For 9 = 45®, however, the discrepancy 
between the ultrarelativistic and mildly relativistic results, shown in Figure 1, 
are significantly larger than those obtained for the circular orbits. For example, 
for 9 = 45®, if v/v c = 2, F ur /F£ 8 (Figure 3, Ramaty, 1968), whereas, ac- 
cording to Figure 1, at the same value of v/v c , j UR y/j % 60. The larger dis- 
crepancy for small values of 9 must be attributed to the combined effects of the 
non- relativistic nature of the particles, the helical orbits and the ionized medium. 
The subsequent calculations, therefore, will entirely be based on the exact 
formalism developed in the previous sections. 

We next consider an ensemble of electrons with an isotropic pitch-angle 
distribution, g (4>) = l/4w , and an energy spectrum of the form 


u(7) ~ (7- l) -r 


i 

I 


(40) 
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with a low-energy cutoff at y = 1.2. It was found (e.g. Holt and Ramaty, 1%9) 
that such a distribution, with r 3, can reproduce the x-ray spectrum between 
100 to 500 keV that was observed for the solar flare of 7 July 1966. At higher 
energies, direct observations of solar flare electrons (Cline and McDonald, 1968) 
also indicate that the spectrum is a power law with an index which may vary from 
flare to flare, but in general has a value close to that deduced from the x-ray 
data. In the subsequent calculations, we shall use Equation (40) withP = 3. 

The emissivities, j ± , and absorption coefficients K ± , as functions of v/ v B , 
are shown in Figure 2, for 0 = 45°, and for c. = 3 and 0.5. As can be seen, for 
the small value of a, the emissivity and absorption coefficients, in both modes, 
are strongly suppressed at low frequencies. For a = 3, this suppression effect 
dissappears, and, at low frequencies, one can clearly see the harmonic structure 
of both j andK. From Figure 2 we also see that for both a - 3 and a = 0.25, 
j . > j + , K_ > K + , but j + /K + > j _/k_ . In an optically thin source, therefore, the 
extraordinary mode will dominate, whereas if the source is optically thick, the 
ordinary mode radiation will be larger than that in the extraordinary mode. 

This is a special result for non- thermal sources, since in an optically thick 
thermal source, because of Kirchhoff's law, j/K will be the same in both modes. 
The dominance of the ordinary mode in a optically thick non-thermal source is 
caused by the greater relative contribution of higher energy electrons to radia- 
tion in this mode. This effect is most pronounced at typical gyro-synchrotron 
frequencies, since in the limiting case of radiation from ultrarelativistic elec- 
trons, the contributions from both modes become equal. 
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As can be seen from Figure 2, the conditions for a given source to be opti- 
cally thin or thick depend critically on the values of the parameters N/BA and a . 

In order to investigate these conditions and the simultaneous dependence of the 
spectrum and polarization on N/BA and a we have evaluated Equations (34) and 
(35), for various values of these parameters, and for 6 - 45°, r = 3, and 
g(<£) = 1/477. The total intensity, I, as a function of v/v B , is shown in Figures 
(3), (4), (5) and (6) for N/BA = 10 9 , 10 u , 10 13 and 10 15 respectively. The 
asterisks on the spectra denote the frequency, v m , at which I + = I_. Since for 

the frequency range of interest (v > v and v > v ) the absolute values of both a , 

p * e 

and a _ are close to unity (a + £ 1; a _ £ -1), the frequency v* also corresponds 

0 0 6 

to complete depolarization (Q = U = V = 0). For v « v * , the source is optically 
thick and I + > I_ . For v » v* the source becomes optically thin, and I_ > I + . Holt 
and Ramaty (1969) made a study of microwave emission from solar flares based 
on the theory developed in the present paper. They showed that the transition 
from the optically thick to the optically thin regime, discussed above could ac- 
count for the observed reversal of polarization of solar microwave bursts. 

From Figures (3), (4), (5) and (6), we can also see that for a given value of 
N/BA one can define a critical value of a, a*, such that: for a < a* the source is 
optically thin, the radiation at low frequencies is dominated by the Razin effect 
and it is polarized in a sense corresponding to the extraordinary mode; for 
a > a* , the radiation at low frequencies is dominated by reabsorption and, at 
v - v* , the polarization changes from the ordinary to the extraordinary mode. 
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In Figure (7) we plot N/BA as a function of a* . This curve was obtained from 
Figures (3), (4), (5) and ( 6 ) as well as from similar calculations for other values 
of N/BA. As indicated, the area to the left of the curve corresponds to an op- 
tically thin source, while for values of N/BA and a corresponding to the area to 
the right of the curve, the polarization reverses at v = v* . 

In order to illustrate the polarization properties of the source, we evaluated 
the Stokes parameters I , Q and V, using Equations (35). (As discussed above, 
because of Faraday rotation, U = 0.) These are shown in Figure 8 as functions 
of i-'A'g for N/BA = 10 13 and for a = 1 and 0.25. As can be seen, for a = 1 , the 
polarization reverses at v/v^z 10. For a = 0,25, this reversal does not occur. 
Since Q « V , the degree of polarization is given approximately by n X V/I . The 
value of 77 is of the order of 10 % to 20 % in the optically thick regime and may 
become as large as 40% at frequencies where the source is optically thin. Since 
Q is very small, the polarization is almost circular. Nevertheless, the measure- 
ment of the sign of Q is possible in principle (see Bekefi (1966) pp. 23). If on#* 
adopts the model of Holt and Ramaty (1969) for the polarization reversal of 
microwave bursts, the measurement of the sign of Q and of the direction in 
which its value is maximal would determine the orientation of the polarization 
ellipse and, therefore, also of B x at the source of the microwave emission. 

We proceed now to consider the effects of an anisotropic pitch-angle dis- 
tribution. For an isotropic distribution, a variation in the angle of observation 
8 has the net effect of reducing the total emission and the emission frequencies, 
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both approximately by a factor of sin 0 . As long as 8 is not equal to or very 
close to zero, these effects are small. On the other hand, for a highly anisotropic 
distribution, the radiation pattern will depend strongly on the angle between the 
direction of observation and that of maximum anisotropy. In order to illustrate 
this effect we have considered a pitch-angle distribution of the form 

g(0) ~ sin m 0 (41) 

and have choeenm = 45 for which g(80°)/g(90®) = 0.5. Using Equations (40) and 
(41), we have evaluated the total intensity, I , for a variety of values of the 
parameters a, 6 and N/BA. From these calculations we find that, as for an 
isotropic distribution, the spectrum at low frequencies will be determined either 
by the Razin effect or by reabsorption, but will essentially be independent of the 
value of i 9 . At the high frequencies, however, the intensity will depend critically 
on the direction of observation. This is illustrated in Figures (9) and (If), for 
a = 1 and 0.25, respectively. 

As in the previous figures, the asterisks denote the frequency at which 
I + = I_ . By comparing the 9 - 90® curve in Figure (9) with that corresponding 
to a = 3 in Figure 5 , we see that if the direction of observation coincides with 
that of maximum anisotropy, the total intensity, spectral shape and frequency 
of polarization reversal for the anisotropic distribution are all comparable to 
the corresponding quantities in the isotropic case (since for the latter 8 t 90®, 
both the intensity and frequency scale are reduced by approximately a factor of 


28 



sin 6). From Figures (9) and (10) we also see that if d <v/ 2, the emission in- 
tensity at high frequencies is strongly suppressed. As discussed above, this 
high-energy cutoff is the direct result of the different emission patterns of 
mildly relativistic and ultrarelativistic electrons. In Figure (10) we see the 
combined effects of the anisotropy and the Razin effect. Since now the radiation 
is suppressed at both the high and low frequencies, a situation may develop in 
which the total emission is also strongly suppressed. 

Finally, we consider the polarization of the spectra shown in Figures (9) 
and (10). By comparing Figure (9) with Figures (3), (4), (5) and "?), we see that, 
for comparable values of a , the frequency of polarization reversal is lower for 
the anisotropic distribution and it decreases with decreasing 8. This is a j ? ’ct 
consequence of the suppression of the emission of the relativistic electrons, 
since as discussed above, the dominance of the ordinary mode in an optically 
thick source results from the enhanced contribution of the high energy electrons 
to the emission at a given frequency. Because of the strong Razin effect, the 
spectra shown in Figure (10) are all due to an optically thin source and their 
polarization corresponds to that of the extraordinary mode. 

Finally, we consider the question of negative absorption. As mentioned 
above, Bekefi, Hirshfield and Brown (1961) have investigated the possibility of 
maser action in a gyro-synchrotron source for which the distribution function f, 
at the appropriate energies, satisfies d f/dy > 0. These authors have considered 
an isotropic electron distribution of the form f ~ (P exp (-b 0 1 ) , where both p 
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and b are non-negative constants, and showed that for p > 0 maser action will 
occur just below v B and its second harmonic. In the present paper we shall in- 
vestigate negative absorption that results from an anisotropic electron distribu- 
tion in a magnetoactive plasma. We shall only consider electron distributions 
for which d i/dy < 0 , so that maser action will result only from the anisotropy. 

We have chosen two specific examples: 

u(y) ~ exp (- 7 ) ( 42 ) 


and 


U(7) ~ y'2.5 (43) 

and a pitch-angle distribution of the form given by Equation (41). Using Equa- 
tions (22) through (25) (for d - n/ 2), we have evaluated K + and K_ for m = 5 and 
for various values of a . We find that for such a moderate anisotropy there is 
no maser action for a - 00 (emission in vacuum). The re sis on for this is that the 
high-energy electrons contribute large positive quantities to both K + and K_ , and 
these cancel the negative contribution that results from the suiisotropy. As a is 
decreased, however, K + becomes negative, in essentially the same frequency 
range (G ) as found by Bekefi, Hirshfield and Brown (1961) for the isotropic 
distribution with an inverted energy spectrum. This is illustrated in Figures 
(11) sind (12) for the spectra given by Equations (42) and (43), respectively. The 
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occurrence of an instability in a given range of frequencies in the presence of a 
plasma, is the result of the suppression of the positive contribution of the rela- 
tivistic electrons to radiation at that frequency. Qualitatively, this result is the 
same as that obtained by McCray (1966) and Zheleznyakov (1967), who showed 
that even though synchrotron radiation in a vacuum can never be amplified, in 
the presence of a plasma maser action is possible. Furthermore, by comparing 
Figures (11) and (12) we also see that the range of frequencies in which K + is 
negative is larger for the exponential spectrum than for the power law. This 
also is the result of the different contributions of the high-energy electrons, 
since the spectrum given by Equation (42), obviously consists of less energetic 
particles than the one given by Equation (43). 

For a moderate anisotropy, such as that given by m = 5, maser action is 
confined to v< v and unless these frequencies can escape from the source, 
amplified radiation will not be observed. As a is further decreased, therefore, 
v will eventually become larger than v B , and, as can be seen from Figures (11) 
and (12), K + will be positive for all frequencies of interest. By the same token, 
since the characteristic frequency , defined above is larger than v , for the 
distributions used in Figures (11) and (12), K_ is positive for all values of v that 
escape from the source. For much larger anisotropies or for inverted energy 
spectra, maser action is not limited to v £ v B , and, therefore, negative values 
of K_ may occur as well. In the present paper, however, we shall not treat 
maser action resulting from such highly anisotropic distribution. 
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Limiting ourself to the spectrum given by Equation (43) and to small values 
of m, we consider the effects of the anisotropy. In Figure (13), K + is plotted for 
a = 2 and for various values of the parameter m. As can be seen, there is no 
maser action for the isotropic distribution but as m is increased, an instability 
sets in and the range of frequencies over which K + is negative increases with 
increasing m. 

Another important aspect of amplified gyro-synchrotron radiation is its 
dependence on the direction of observation. In order to investigate this effect, 
we have evaluatdff the functions G and H for various values of the observation 
angle 6. We have used Equations (18) through (21) for 6 / -n/ 2 and Equations 
(22) through (25) for 6 - n/ 2 . The total intensity I , obtained from Equations 
(34) and (35) is plotted in Figure (14) for several values of 9. As can be seen, 
the radiation is highly amplified at v < , but falls off rapidly as the angle be- 

tween the direction of observation and that of maximum anisotropy is increased. 
Since the amplified radiation is entirely due to the ordinary mode and since it is 
confined t o9Z n/ 2, even though Faraday rotation is large, the emission is 100% 
linearly polarized. Because of its highly collimated nature and large degree of 
polarization, this type of radiation may be responsible for radio emission from 
the recently observed pulsars. 

V. SUMMARY 

We have developed a detailed formalism for the calculation of the emissivity 
and absorption coefficient, for a given direction of observation, from an ensemble 
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of electrons with arbitrary energy and pitch-angle distributions in a cold col- 
lisionless electron plasma and a static uniform magnetic field. By assuming 
that all spatial distributions are homogeneous, we have solved the transfer 
equations in the limit of large Faraday rotation, i.e. when the phase difference 
between the ordinary and extraordinary modes is completely random so that 
the two modes propagate independently. Tula condition is generally well satis- 
fied for gyro- synchrotron sources. 

Using this formalism we have studied the influence of reabsorption anu the 
Razin effect on the spectrum and polar zation of the radiation. Both effects tend 
to lower the flux at low frequencies, but since the Razin effect suppresses not 
only the emissivity but also the absorption coefficient, a source which is opti- 
cally thick in the absence of an ambient plasma, may become optically thin if 
the influence of the medium is taken into account. Gyro-synchrotron radiation 
from an optically thin source is polarized in the extraordinary mode, but be- 
cause of the non-thermal nature of the radiating electrons, in an optically thick 
source the ordinary mode dominates. This transition is characteristic of gyro- 
synchrotron radiation and it may account for the observed polarization reversal 
of solar microwave bursts. 

We have also investigated the effects of an anisotropic electron distribution. 
For such distributions, if the direction of observation differs significantly from 
that of maximum anisotropy, the radiation at high frequencies is strongly sup- 
pressed. This cutoff, together with the low frequency Razin suppression can 
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effectively suppress the entire radiation spectrum. Moreover, the pitch-angle 
anisotropy, by suppressing the radiation from the high-energy electrons, has 
the additional effect of lowering the frequency of polarization reversal, mentioned 
above. 

Finally, we have studied the possibility of maser action. We find that an 
anisotropic distribution can cause negative absorption, even if the electron 
-spectrum, at all energies, decreases with increasing energy. For moderate 
anisotropies maser action will be limited to the ordinary mode and to fre- 
quencies just below the gyro frequency. The radiation will be 100% polarized 
and if the pitch- angle distribution is peaked around circular orbits, the ampli- 
fied radiation will be linearly polarized and confined to a narrow cone centered 
around the direction of maximum anisotropy. 
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FIGURE CAPTIONS 


Figure 1. 


Figure 2. 


Figure 3 . 


Figure 4. 


Figure 5. 


Figure 6. 


Figure 7. 


Gyro- synchrotron emission spectra from monoenergetic electrons. 
The dashed lines were obtained from the ultrarelativistic approxima- 
tion (Equation 37) and the solid lines from the exact formalism. 

The emissivity (solid line) and absorption coefficient (dashed line) 
from an ensemble of electrons with isotropic pitch-angle distribution 
and power law energy spectrum. 

The radiation intensity from an isotropic, power law electron dis- 
tribution for various values of the Razin parameter a and the re- 
absorption parameter N/BA equal to 10 9 . 

The radiation intensity from an isotropic, power law electron dis- 
tribution for various values of the Razin parameter a and the re- 
absorption parameter N/BA equal to 10 u . 

The radiation intensity from an isotropic, power law electron dis- 
tribution for various values of the Razin parameter a and the re- 
absorption parameter N/BA equal to 10 13 . 

The radiation intensity from an isotropic, power law electron dis- 
tribution for various values of the Razin parameter a and the re- 
absorption parameter N/BA equal to 10 15 . 

Relationship between the Razin and reabsorption parameters, a and 
N/BA, respectively. The area to the left corresponds to optically 


Figure 8. 


Figure 9. 


Figure 10. 


Figure 11. 


Figure 12. 


Figure 13. 


Figure 14. 


thin sources, whereas in the area to the right, the source is optically 
thick below v* but becomes optically thin above v* . 

Variation of the Stokes parameters I , Q and v with frequency for 
two choices of a and N/BA = 10 13 . 

The radiation intensity from a highly anisotropic electron distribu- 
tion for various values of the emission angle 6 and for a = 1. 

The radiation intensity from a highly anisotropic electron distribu- 
tion for various values of the emission angle 6 and for a = 0.25. 
Absorption coefficients in the ordinary mode for various values of 
the Razin parameter a and a power law electron spectrum. 
Absorption coefficients in the ordinary mode for various values of 
the Razin parameter a and an exponential electron spectrum. 
Absorption coefficients in the ordinary mode for various values of 
the anisotropy parameter m . 

Total radiation intensity which exhibits a strong enhancement due to 
maser action at v < v B . 
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